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Abstract 

The Zipoy-Voorhees family of static, axisymmetric vacuum solutions 
forms an interesting family in that it contains the Schwarzschild black 
hole excepting which all other members have naked singularity. We ana- 
lyze some properties of the region near singularity by studying a natural 
family of 2-surfaces. We establish that these have the topology of the 
2-sphere by an application of the Gauss-Bonnet theorem. By computing 
the area, we establish that the singular region is 'point-like'. Isometric 
embedding of these surfaces in the three dimensional Euclidean space is 
used to distinguish the two types of deviations from spherical symmetry. 
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I. INTRODUCTION 



The Weyl class of static, axisymmetric vacuum solutions of Einstein equation is very 
well studied and a systematic procedure for constructing solutions is available [jXJ] . This 
class of solutions can be looked upon as describing finite size static bodies with axial 
symmetry. These admit a regular horizon only when the symmetry turns spherical and 
the solution is the Schwarzschild black hole. Among these is a class of so called 'prolate 
spheroidal' solutions found by Zipoy 0. A special family of this class is the subject of 
this note. 

Solutions of this family are parameterized by two parameters, m and A. All are 
asymptotically flat and all except A = 1, which is the Schwarzschild solution, have 
'naked singularity'. It is interesting to note that arbitrarily small deviations from the 
Schwarzschild black hole value of A = 1, which also introduces deviations from spheric- 
ity, convert the event horizon into an eternal naked (curvature) singularity. This makes 
the study of this family interesting and this is undertaken here. 

In section II, we sketch an alternative derivation of the solution. Traditional ap- 
proaches to construction of Weyl solutions have involved use of 'Newtonian potentials' 
corresponding to various source configurations eg [|],[|. We discovered this solution while 
studying the Kerr metric derivation given by Chandrasekhar |4]] in some other context. 
The solution is arrived at essentially by specializing Chandrasekhar's equations 0] lead- 
ing to the Kerr solution, by switching off rotation (u = 0). This is the static limit of 
the Kerr geometry which would reduce to the Schwarzschild black hole by the demand 
of the existence of a regular horizon. 

In section III, we analyze some properties of the solution. After stating some ele- 
mentary properties we focus on its t =constant, r =constant surfaces, S 2 . We analyze 
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its topology and intrinsic geometry and examine its 'shape' as viewed in a three dimen- 
sional Euclidean space. It turns out that for the solutions for A > 1 and r very close to 
2m, these surfaces are not completely embeddable in the Euclidean three dimensional 
space. This is a rather unusual and interesting feature of this family of solutions. What 
physical implication this feature has is at present an open question. 



In the final section, we comment on possible physical implications of the results and 
open issues. The notation and conventions used are those of Chandrasekhar |4[] with the 
metric signature (H ). 



II. DERIVATION OF THE SOLUTION 



The solution is fairly straightforward to derive. We begin by the general equations 
given in M for stationary, axisymmetric space-times and specialize to the static form of 
the ansatz by setting u = 0. The organization of equations and subsequent simplifica- 
tions is identical to that given in [|J. Hence we give only the key steps. 



The metric is given by, 

The metric coefficient functions depend only on r, 6 and A =exp[2(/i3 — z^)] is freely 
specifiable. Defining /3 = ip + v, the vacuum equations become, 

d r (VAd r eP) + d e {-j=d e eP) (2) 



d r {e^d r {4, - i/)} + d {-=d e {il> - i/)} (3) 

v A 



d r d e p - d r f3d e {i 2 - d d (3d r [i 3 + d r vd Q v + d r ^d e ^ (4) 
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I— 4 
4VA{<9 r /?<9 r /i 3 + d r i)d r v} ={0 e /?%i2 + do4>d e v} 

v A 



2e- (3 {d r (^d r e f3 ) - d e (^=d e e^} (5) 

v A 



In order to admit possibility of a horizon, following Chandrasekhar, we set 



exp{/3(r,6»)} = JA(r)f(9). Equation (2) then leads to A = r 2 + d x r + d 2 , f{6) = sin(9). 
In anticipation of spherical topology we have required f(9) to vanish at the 'poles' at 
9 = 0, 7r. We will assume that di constants are such that A = could have two real 
roots and we may write A(r) = (r — r + )(r — rS). The remaining equations now become 
equations for v and \xi (say) by eliminating ip in favour of u, j3 and by eliminating fi 3 in 
favour of /X2, A. The class of solutions we are interested in is obtained by assuming that 
u(r, 9) is independent of 9. The equations simplify further and lead to the following 
solution: 



C n r r-r +1 2C n ( sin 2 9 4 
in { ±} _ N2^{— — + ^. 



/i 3 = // 2 + -£nA , A = (r - r + )(r - r_) 

^fo(sin(^)) + ^n(A) , if) = p - u (6) 

Here C, £) and $o (and of course r±) are constants of integration. By redefining 
r — > r + r_ one sees that r_ now appears only in the combination r + — r_ = 2m > 0. For 
vacuum solution, the metric is always defined up to an over all multiplicative constant. 
We also have the freedom to redefine t and 0. Using these one can effectively absorb the 
constants D, $o- Trading C for mA, the final form of the metric can be written as: 

ds 2 = F A dt 2 - F~ A G~ B dr 2 - AF~ A G~ B d9 2 - AF- A sin 2 (9)d^ 2 (7) 

where, A = r 2 - 2mr , F = — , G = 1 + —sin 2 (9) and 5 = A 2 - 1 

r A 

Let us quickly remark that the solution can be expressed in a standard form [|TJ by 
the following identifications: 



p = y/AstnO , z = -^cosO , e 2U = F A and e 2k = G~ A \ (8) 

2 dr 

So that the metric is given by, 

ds i = e 2U dt 2 -p 2 e- 2U d<p 2 -e- 2U+2k (dp 2 + dz 2 ) (9) 

We also note that this solution corresponds to the 'prolate' family of solutions found 
by Zipoy @. In equations (43-45) of j|, put 7 = 0. A trivial shift of his r gives the 
above solution with his (3 same as our A. Voorhees || also obtained the solution as an 
illustration of his method of constructing Weyl solutions. He also elaborated some of its 
properties. 

In the next section we recall some of the basic properties and then focus on the t = 
constant and r = constant surfaces S 2 . 

III. PROPERTIES OF THE SOLUTION 

We have a two parameter family of static, axisymmetric solutions potentially con- 
taining a horizon. For A = 1 we obtain the Schwarzschild solution with mass m, r 
being the usual areal radial coordinate and r = 2m being the horizon. The spherical 
symmetry is restored for this value. For large r, the behaviour of the solution shows that 
it is asymptotically flat with mA being the ADM mass. We are primarily interested in 
parameter values closer to the Schwarzschild case. Thus we take m > and A > so 
that the mass is also positive (m = gives Minkowski space-time). We now concentrate 
on A / 1. 

To locate possible singular regions, we compute curvature components. For example 
(in Chandrasekhar's notation), 

R 0W1 = -AG B A A - 2 r- 2A mr(l - ^ A + 1 ) m ) (10) 



5 



As r — > 2m, A — > 0. For A 7^ 1 (B ^ 0), the net power of A becomes A — A 2 — I which 
is always negative for all rea/ A. Thus we get a curvature singularity as A — > 0. Other 
curvature components and Kretschmann scalar show similar behaviour. For A = 1, 
there is no contribution from G and the curvature component reduces to — ^. Since for 
r — > 2m itself we encounter curvature singularity, there is no point in exploring r < 2m 
region. Since there is no other r > 2m where the norm of the stationary Killing vector 
(gtt) can vanish, the singularity must be naked. Thus except for A = 1 where one has a 
non singular Killing horizon, for all other values of A the space-time has naked singular- 
ity (diverging curvatures). This is an eternal naked singularity in an asymptotically flat 
space-time with positive ADM mass. From now on we limit ourselves to the non-singular 
coordinate range r > 2m. In order to gain some intuitive understanding of the region 
near the singularity as well as to see how the A ^ 'effects' are seen in large r regions, 
we focus now on the t = constant and r = constant surfaces, generically denoted by S 2 . 
This naturally leads us to study the intrinsic geometry of E 2 . 

The metric as given has all the coordinates as local coordinates. We have implicitly 
assumed the 9, <fi to be the standard spherical polar angles on S 2 . As such the non- 
singular form of the metric is not directly in conflict with such an interpretation of the 
coordinates. (While can be taken as azimuthal angle referring to the axisymmetry, in- 
terpretation of 9 as the polar angle is not automatic.) One way to check that £ 2 indeed 
can have spherical topology is to verify the Gauss-Bonnet theorem. This is precisely 
what is done below. 

It is straight forward to compute the Ricci scalar for the intrinsic metric on a £ 2 . 
Since r is a constant, defining a = AF~ A and ft = ^ one obtains the Ricci scalar as, 




(11) 



-{1 + (3sin 2 9Y B - l \{(5B - 1) - (3{1 + B)sin 2 6} 



(12) 



a 
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If topology of £ 2 is indeed S 2 , then integral of \Jdet(g)R(g) over the full range of 
the coordinates should equal 871. In terms of x — cos9, the integral can be expressed as, 

, 1 1+/3 _ (-Q _|_ ]_W2 

- 1 = ^ & ( ^"^ 2)1 _ B/2 (13) 

Clearly for 5 = (Schwarzschild), the right hand side is 1 as expected. Defining 
£=(l + /3)//3, we get, 

i- ^(i)R(g) = (£ - 1)- B / 2 [(S + 1) - 2^ 5 ] <?(£, 5) (14) 

where, <?(£,£)= f dx - x 2 ) B ' 2 

Jo 

Note that £ > 1 and hence the function #(£, _B) is well defined and is in fact differ- 
entiable. One can simplify the £ dependence by defining x = \/^y in the integral. It 
is now easy to evaluate the ^-derivative and check explicitly that the right hand side is 
indeed 1 for all values of B and (3. Thus topology o/E 2 surfaces is indeed that of S 2 . 

Next we compute the areas of these spheres. This will allow us to estimate the 'size' 
of the region of high curvature. The area is given by, 



Area(r) = 4vr r 2A A (1 ~ A)2/2 [ [' dx {A + m 2 (l - x 2 )}^ 

Jo 



-A 2 )/2 



(15) 



Putting r = m(/i + 1), x — fiy and 7 = (1 — A 2 )/2, the area can be expressed as, 

Area(/i, A) = 4tt m 2 ^ (fx + 1)( 1 + A ) 2 / 2 (^ - 1)(^) 2 / 2 f 1/fl dy (1 _ y 2y (16) 

Jo 

Since r > 2m we have /i > 1 and manifestly 7 < 1/2. The area is clearly well defined 
and is a positive, finite quantity for all finite /i > 1. Our interest is to estimate the 
behaviour of the area as /i — > 1. Once again for A = 1, we recover the expected answer. 

It is easy to see that for < A < 2, A ^ 1, the area vanishes as r — > 2m while for 
larger values of A the area blows up. The vanishing area indicates that the singular re- 
gion is 'point-like'. For diverging area we can not say so. Since we are mainly interested 



in values of A near the Schwarzschild value, we see that for these values the area vanishes 
and hence the curvature singularity is 'point-like'. The coordinate value r = 2m really 
corresponds to areal radial coordinate vanishing. For the subsequent analysis we will 
limit ourselves to solutions with 'point-like' singular region, i.e. to < A < 2. 

Note: If we assume that the Ricci scalar satisfies a < R < b everywhere on a spherical 
surface for some constants a, b depending on the surface, then the Gauss-Bonnet integral 
for spherical topology gives, 

a (Area) < 8tt < b (Area) (17) 

Thus if the area vanishes then the Ricci scalar must be unbounded above while if the 
area diverges then either the Ricci scalar must vanish in a precise manner (eg usual S 2 
metric for large radius) or a must be negative or a can go to — oo. One can check from 
our expressions that this is indeed the case. Observe that, the A — > behaviour of the 
Ricci scalar is different along the poles and along generic 9 directions. 

Having estimated the 'size' of the singular region, we now look for its 'shape'. All 
our intuitive understanding of 'shape' of an object is based on its embedding in three 
dimensional Euclidean space. A natural formulation of 'shape' of S 2 is to look for its 
embedding in the Euclidean space of dimension 3. The embedding is to be such that 
the induced metric on the image is the same as the intrinsic metric on it, in short, an 
isometric embedding. 

Using the symmetries of the Euclidean space and choosing the axis of symmetry to 
be the Z-axis, a natural ansatz for the embedding is: 

X(e,4>) = x{9) cos(f) , Y{9,<j>) = y{9) sin<t> , Z(9,4>) = z{9) . (18) 

The demand that the induced metric be axisymmetric (^-independent) and diagonal, 
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leads to y{9) = ±x(9) and without loss of generality we choose the + sign. Equating 
with the intrinsic metric then gives (prime denotes derivative with respect 9), 

x{9) = ±v^ sin6 , (z') = ± y/a [(1 + (3 sin 2 6)- B - cos 2 6] 1/2 (19) 

Without loss of generality we can choose positive sign for x{6) and negative sign for 
z{9) equations above. The equation for z{9) is an ordinary differential equation with one 
constant of integration which corresponds to the choice of origin. It is invariant under 
reflection about the equator. We can thus solve the equation for < 6 < n/2. For 
convenience, let us denote the expression in the square brackets in the above equation 
by f(9,P,B). We are interested in the behaviour of this function for < 9 < ir/2, for 
< (3 and for -1 < B < 3 (0 < A < 2). 



Clearly / is non-negative for all (3, 9 if B < 0. however for B > and for some (3, it 
can be negative for 9 near the poles, indicating failure of complete embedding. 



Computing the first derivative of / with respect to 9 one sees that 9 = 0,7r/2 are 
two extrema with a possibility of a third extremum at some 9. This third extremum, if 
exists, is always a minimum. For j3B < 1, the third extremum does not exist while for 
fiB = 1, it coincides with 9 = 0. Then / is non-negative for all 9. Thus for f}B < 1, 
complete embedding is possible. For f3B > 1, the minimum at 9 > exists and then / 
is negative for < 9 < 9 max . The 9 max is to be obtained by solving the transcendental 
equation / = 0. Thus for (3B > 1, caps nears the poles are not embeddable. 



The condition (3B > 1 corresponds to r near 2m. Since for B close to zero, (3 is very 
large and thus corresponds to r very close to 2m. The farthest r would be for the largest 
B (= A 2 — 1) which for us is 3 (A < 2). This translates to the farthest r being 3m. 
Thus for r > 3m, we will always have complete embedding. Curiously, r = 3m happens 
to be the radius of the photon circular orbit for the Schwarzschild black hole. 
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While precise shape for complete embedding can be plotted by numerically integrat- 
ing the equation, that is however not our main concern. In general the 'shape' will be 
an ellipsoid. For B close to zero, the ellipsoid will shrink as A — > 0. 

IV. DISCUSSION 

One of the striking features of this family of solutions is that for a slightest devia- 
tion from spherical symmetry, the horizon disappears and instead a naked singularity 
appears which is eternal by virtue of the staticity of the space-time. The departure from 
sphericity, equivalently from the Schwarzschild geometry, is indicated by the parameter 
A ^ 1. Note that A = +1 corresponds to the Schwarzschild black hole while A = — 1 
corresponds to the negative mass Schwarzschild solution. We are assuming of course 
that the parameter m is positive. In order to include the Schwarzschild black hole, we 
have chosen A > 0. The parameter value A = 2, demarcates whether the area of E 2 for r 
near 2m vanishes (A < 2) or diverges (A > 2). Likewise, the Schwarzschild value A = 1 
demarcates whether the surface £2 for r near 2m is completely embeddable (A < 1) or 
is partially embeddable (A > 1). 

Of course it is more likely that the solutions are valid for r sufficiently large and that 
a physical body will fill in the interior region. This leads one to the question of a source 
configuration to which the above family is an exterior space-time . While it would be 
natural to expect that a two surface across which an interior solution is matched to the 
above solution is indeed one of the £2 surfaces, it is by no means obvious that it must 
be so. For instance one could have a 2-surfaces defined by r = r{9). The Killing norm 
however will not be constant over such surfaces. For large r values this will correspond 
to the matching surface not having constant Newtonian potential. One could of course 
analyze the shape of any such matching surfaces by similar techniques. One must em- 
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phasize that the Euclidean space used for the isometric embedding is not the physical 
space whose metric near the vicinity of the surface is not Euclidean. It has been used 
as a way to discriminate among different members of the family []. 

These solutions could be viewed as a subset (static) of possible end states of a non- 
spherical collapse. Black hole uniqueness results would indicate that if censorship holds 
then either spherical symmetry must be restored (staticity results) or rotation must set 
in (only stationarity results) during the collapse. For the possibility of static end state, 
it is only a single special value for which a naked singularity is avoided. In the light of 
these observations, it would be exceedingly interesting to see what features emerge in 
an 'axisymmetric' collapse with the above family being used for exterior matching just 
as the Schwarzschild is used in the study of spherical collapse. 
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1 While the paper was being written, a preprint by N. Pelavas, N. Neary and K. Lake appeared 
on the e-print archive which applies similar methods to study properties of the instantaneous 



ergo surface of the Kerr metric, gr-qc/0012052. 
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